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Spectroscopy of geo-neutrinos from 2056 days of Borexino data 
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We report an improved geo-neutrino measurement with Borexino from 2056 days of data taking. 

The present exposure is (5.5 ± 0.3) x 10®® protonxyr. Assuming a chondritic Th/U mass ratio of 
3.9, we obtain 23.7 (i^'^istat) /o'6('52/s) geo-neutrino events. The null observation of geo-neutrinos 
with Borexino alone has a probability of 3.6 x 10“® (5.9 (j). A geo-neutrino signal from the mantle is 
obtained at 98% C.L. The radiogenic heat production for U and Th from the present best-fit result 
is restricted to the range 23-36 TW, taking into account the uncertainty on the distribution of heat 
producing elements inside the Earth. 

PACS numbers: 13.35.Hb, 14.60.St, 26.65.-l-t, 95.55.Vj, 29.40.Me 


Geo-neutrinos are electron anti-neutrinos {Oe) pro¬ 
duced by (3 decays of long-lived isotopes, which are nat¬ 
urally present in the interior of the Earth, such as decays 
in the chains, and ®®K [TJ [2]. Geo-Pe 


measurements have been reported by Borexino nil] and 
KamLAND [S1E|- Here we present improved reactor neu¬ 
trino and geo-Pe measurements performed by Borexino. 
Borexino is an unsegmented liquid scintillator detector in 
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operation at the underground Gran Sasso National Lab¬ 
oratory, Italy [ 8 ]. In liquid scintillator detectors be are 
detected via the inverse (3 decay process, Ve+p e^-l-n, 
with a threshold of 1.806 MeV. The deposited be energy 
results in a prompt signal induced by the positron, and 
includes annihilation photons. The visible energy is re¬ 
lated to the be energy as Eyis = — 0.784 MeV. A 

delayed signal induced by the neutron capture on pro¬ 
tons produces the 2.22 MeV gamma-ray, providing a de¬ 
layed coincidence signal with a mean capture time of 
259.7±1.3(stat)±2.0(syst) ps [7]. In Borexino 278 tons 
of ultra-pure organic liquid scintillator (pseudocumene 
doped with 1.5 g/1 of diphenyloxazole) is confined within 
a thin, spherical nylon inner vessel (IV) with a nominal 
radius of 4.25 m and concentrically placed inside a stain¬ 
less steel sphere (SSS) with radius of 13.7 m and equipped 
with 2212 8 ” photomultipliers (PMTs). The volume be¬ 
tween the SSS and the IV, divided in two regions by a 
second outer nylon vessel (OV) is filled with 890 tons of 
pseudocumene with 3 g/1 of light quencher dimethylph- 
thalate to shield the core of the detector against 7 radi¬ 
ation and radon from the PMTs and the SSS. The SSS 
is immersed in an water tank (WT) instrumented with 
208 PMTs as a Cerenkov active muon veto, also acting 
as an shield against 7 -rays and and neutrons from the 
surrounding rock. 

In Borexino, known sources of be events are nuclear 
reactors and geo-neutrinos. For geo-neutrinos, 

^^^Th are the only isotopes abundant enough to signifi¬ 
cantly contribute events in Borexino above the inverse /3 
decay threshold. The light yield in Borexino is measured 
to be about 500 photoelectrons (p.e.) / MeV and the en¬ 
ergy resolution scales as ~5%/-\/(i?) [S]. The Borexino 
liquid scintillator shows a high pulse shape discrimination 
efficiency in separating high ionizing particles (protons, 
a’s) from electrons and gamma-rays [9j. 

The data reported here were collected between Decem¬ 
ber 15, 2007 and March 8 , 2015 for a total of 2055.9 days 
before any selection cut. The number of PMTs for the 
present data set has been declining with time, from 1931 
to 1525 (average 1730), with small run by run fluctua¬ 
tions. As reported in [3] the event energy is a calibrated 
non-linear function of the number of detected p.e. We 
perform the analysis directly in number of p.e. We dis¬ 
card events occurring within 2 ms of every muon crossing 
the outer detector and within 2 s of muons crossing the 
inner detector to reject neutrons and long-lived cosmo- 
genic radioactivity, respectively. This cut reduces the 
live-time to 1841.9 days. We apply the following addi¬ 
tional selection cuts: ( 1 ) prompt scintillation light: Qp > 
408 p.e. (i.e. 1.022 MeV corrected for the energy resolu¬ 
tion); ( 2 ) delayed signal scintillation light: 860 < Qd < 
1300 p.e. (neutron capture peak) ; (3) correlation dis¬ 
tance between prompt and delayed signals: Ai? < Im; 
(4) correlated time between prompt and delayed signals: 
20 < At < 1280 pLS] (5) pulse shape discrimination with 


Gatti filter uni: 5a/3 < 0.015 for delayed signals; ( 6 ) 
multiplicity cut: selected event are neither preceded or 
followed by neutron-like events within a 2 ms window; 
(7) dynamical fiducial volume [9]: every prompt signal 
has a reconstructed vertex >30 cm away from the time- 
varying IV surface; ( 8 ) FADG cut: independent check 
of candidate events features by a 400 MHz digitizer ac¬ 
quisition system. The combined efficiency of the cuts is 
determined by Monte Garlo to be (84.2±1.5)%. The to¬ 
tal efficiency-corrected exposure for the present data set 
is 907±44 tonxyr. We have identified 77 be candidates 
passing all the selection cuts. 

TABLE I. Estimated backgrounds for be given in number of 
events. Upper limits are given for 90% C.L. 


‘'Li-^He 

0 iq4+U-1^5 

u.iy^_0 089 

Accidental coincidences 

0.221T0.004 

Time correlated 

D r)Qc:+0.029 

U.UOO_Q Q23 

(a,n) in scintillator 

0.165±0.010 

(Q,n) in buffer 

<0.51 

Fast n’s (/r in WT) 

<0.01 

Fast n’s (p in rock) 

< 0.43 

untagged muons 

0.12±0.01 

Fission in PMTs 

0.032±0.003 

214Bi-214po 

0.009T0.013 

Total 

o.78ro:iS 

< 0.65(combined) 


The probability of be-mimicking background events 
leaking into the dataset was evaluated as follows: ( 1 ) 
The rate of accidental coincidences has been searched for 
by shifting the delayed time window to 2-20 seconds and 
keeping all other cuts unchanged. The energy spectrum 
of these events is limited to <3 MeV. (2) Time correlated 
events have been searched for in the (2 ms, 2 s) time win¬ 
dow. A negligible amount of correlated events with a ~1 
s time constant were identified and their contribution in 
the be time window determined. (3) (Q;,n) background for 
be search has been extensively discussed elsewhere Hi]. 
The average rate of ^^°Po in the dataset is determined 
to be (14.1±0.2) counts/(day-ton). (4) Cosmogenic ra¬ 
dioactive isotopes which decay via /3-|-n, namely ®Li-®He 
[7|, have been studied in the (2 ms, 2 s) time window after 
a muon crossing the inner detector. These events have 
a wide energy spectrum with maximum at ~5 MeV. (5) 
^^^Rn in the liquid scintillator can produce background 
events through time-correlated + decays of ^^^Bi 

and This decay sequence has a time constant 

close to the neutron capture time following a be inter¬ 
action. This background has been estimated by individ¬ 
ually tagging ^^^Bi-^^^Po decays. ( 6 ) Backgrounds from 
fast neutrons, untagged muons and spontaneous fission 
decays in the PMTs are the same as in previous papers 
H Hj. Table E summarizes the estimated backgrounds 
for be candidates, expressed in number of events. The 
combined upper limit is obtained by Monte Garlo. The 
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signal-to-background ratio is ~100. 

We have performed an un-binned likelihood fit of the 
energy spectrum of selected prompt candidate events 
0, shown in Fig. The reactor and geo-neutrinos spec¬ 
tra are obtained by Monte Carlo and the backgrounds 
considered in this analysis are reported in Table |l] The 
Monte Carlo spectra have been determined as reported 
in |3]. The reactor neutrinos signal has been calculated 
adopting the data from IAEA m updated to 2014 and 
the method described in [12]. For the first quarter of 2015 
we have used the values from 2014. For the present expo¬ 
sure we predict (87 ± 4) TNU events from nuclear reac¬ 
tors, where 1 Terrestrial Neutrino Unit (TNU) = 1 event 
/ year / 10^^ protons. The log-likelihood function has two 
signal components, Sgeo and Sreact, left free, and three 
background components, SuHe, San, Sacc, constrained 
to the values and errors reported in Tab. |T| These com¬ 
ponents account for 75% of the total background. The 
other components were left out due to the uncertainty in 
their energy spectrum. Combined, they contribute ~1% 
to the best fit and their contribution to the systematic 
uncertainty is absorbed in the uncertainty on the energy 
scale. 

Using the value ratio for the masses of Th and 
U, m(Th)/m{U) = 3.9, suggested by the chondritic 
model, our best fit yields S'geo = 

events (43.5t(|!)'4(sto<)+2'4(s?/s) TNU) and Sreact = 
52.7^7Y(staf)to;g(sys) events (96.6t4|;2('®^a^)^5’o(a2/'®) 
TNU). When expressing the results in TNU, system¬ 
atic uncertainties from both the exposure (4.8%) and the 
Monte Carlo energy calibration (1%) are included. Only 
the Monte Carlo calibration uncertainty is relevant when 
using the number of decays. 

In Fig. [^we show the 1, 3 and 5a contours from the 
log-likelihood fit. Borexino alone observes geo-neutrinos 
with 5.9cr significance (Fig. [^. The null hypothesis for 
geo-neutrino observation has a probability equal to 3.6 x 
10“®. The measured geo-neutrino signal corresponds to 
Pe fluxes at the detector from decays in the U and Th 
chains of (l>{U) = (2.7 ±0.7) x 10® cm“^s“^ and (l>{Th) = 
(2.3 ± 0.6) X 10® cm“^s“^, respectively. Statistical and 
systematic uncertainties are added in quadrature. 

Fig. shows the probability contours obtained by per¬ 
forming the fit leaving the U and Th spectral contribu¬ 
tions as free parameters. The U and Th best-fit contri¬ 
butions are shown in Fig. This measurement shows 
how Borexino, with larger exposure, could separate the 
contributions from U and Th, and demonstrates the abil¬ 
ity of this detection technique to perform real-time spec¬ 
troscopy of geo-neutrinos. 

The radiogenic heat production for U and Th, H{U + 
Th), from the present best-fit result is restricted in the 
range 23-36 TW (see Fig. |^. The range of values in¬ 
cludes the uncertainty on the distribution of heat produc¬ 
ing elements inside the Earth. The model-independent 



Prompt Event Energy [p.e.] 

FIG. 1. Prompt light yield spectrum, in units of photoelec¬ 
trons (p.e.), of i>e candidates and the best-fit. The best-fit 
shows the geo-neutrino and reactor neutrino spectra (dotted 
lines) assuming the chondritic ratio. Colored areas show the 
result of a separate fit with U (blue) and Th (light blue) set 
as free and independent parameters. 

analysis yields a radiogenic heat interval 11-52 TW 
(69% C.L.) for H{U + Th). Adopting the chondritic 
mass ratio above and a potassium-to-uranium mass ra¬ 
tio m{K)/m{U) = 10"^, the total measured terrestrial 
radiogenic power is P{U + Th + K) = 33^20 TW, to 
be compared with the global terrestrial power output 
Ptot =47±2 TW [Hj. 

The contribution to the total geo-neutrino signal from 
the local crust (LOG) is estimated to be S'geo(LOG) = 
(9.7±1.3) TNU [T5|. Gonsidering the contribution from 
the rest of the crust (ROG) [16], the signal from the 
crust in Borexino is calculated as ^geo(LOG±ROG) = 
(23.4±2.8) TNU. In order to estimate the significance of 
a positive signal from the mantle we have determined the 
likelihood of SgeoiMantle) = Sgeo — S'geo(LOG±ROG) 
using the experimental likelihood profile of Sgeo and a 
gaussian approximation for the crust contribution. The 
non-physical region, Sgeo{Mantle) < 0, is excluded. 
This approach gives a signal from the mantle equal to 
SgeoiMantle) = 20.9t4g;3 TNU, with the null hypothe¬ 
sis rejected at 98% G.L.. 

An updated measurement of Ve^s with Borexino is pre¬ 
sented. We show that Borexino-only data measure geo¬ 
neutrinos with 5.9 ct significance. We also shows that 
the background level in Borexino allows to perform a 
real time spectroscopy of geo-neutrinos, currently limited 
only by the size of the detector. 
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(UMO-2012/06/M/ST2/00426). We acknowledge the fi¬ 
nancial support from the UnivEarthS Labex program of 
Sorbonne Paris Git (ANR-lO-LABX-0023 and ANR-11- 
ID EX-0005-02). We acknowledge the generous support 
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5react[TNU] 

FIG. 2. Best-fit contours for 1, 3 and 5a for the statistics 
reported in this paper. 



5£/[TNU] 

FIG. 3. Best-fit contours for 1, 2 and 3 ct for the statistics re¬ 
ported in this paper and for an unbinned likelihood fit with U 
and Th kept as distinct and free parameters. All other param¬ 
eters in the fit were kept unchanged. Dashed line corresponds 
to the chondritic assumption. 

Sasso (LNGS). We acknowledge the collaboration of J. 
Mandula for providing detailed data on reactor load fac¬ 
tors. 
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